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We formulate eikonal approximation to the calculation of high energy scattering amplitude in the 
frame where both colliding objects are very energetic. We express the eikonal scattering matrix 
in terms of the color charge densities of the colliding objects. The calculation is performed in the 
Hamiltonian formalism. We also show that the appearance of the longitudinal electric and magnetic 
fields immediately following the collision is fully taken into account in the eikonal approximation. 

INTRODUCTION. 

Partonic level eikonal approximation has been extensively used in recent years in the study of high energy evolution 
of hadronic scattering [1] . As explained in [2] , this approximation assumes only integrity of partons throughout the 
scattering event rather than the integrity of the whole hadronic state. It thus allows for inelastic rather than only 
elastic final states. The approximation is most convenient to use in the situation where one of the colliding hadrons 
carries almost all the energy of the collision. In this case the picture of the scattering is that of the partons from the 
wave function of the fast hadron propagating through the color fields associated with the target - an almost static 
hadron[3]. This setup is very convenient when the scattering is strongly asymmetric, namely when the "projectile" 
hadron is small and contains a small number of partons, while the target hadron is large and is characterized by 
strong fields. Hamiltonian description in this case is given straightforwardly in terms of the light cone wave function 
of the projectile [3], [4]. 

On the other hand in a more symmetric situation, such as scattering of two small objects or two large objects, like 
nuclei, it is advantageous to describe the scattering event in a symmetric frame where both colliding objects have 
comparable energies. This requires the use of a symmetric gauge condition. In Lagrangian formulation this has been 
done for both "p-p" and "A- A" scattering [5], [6], [7]. We are however not aware of a Hamiltonian formulation of this 
problem in a symmetric frame and gauge. It is desirable to have such a formulation, since Hamiltonian approach 
proved to be very fruitful and intuitive in the study of high energy scattering [4], [8], [9], [10]. 

Our purpose in the first part of this paper is to provide such a Hamiltonian formulation. We use a symmetric gauge 
A$ = with the subsidiary gauge fixing condition A^(t = 0) = and show how to calculate the eikonal S'-matrix in 
the "classical approximation", namely in the situation when the filds of both incoming hadrons are parametrically 
large. 

The second part of this paper is a short comment on the nature of the longitudinal fields that are created in the 
initial stages after the collision. The fact of the appearance of such fields has been realized only in recent years and 
this realization prompted a lot of discussion [11]. It has been proposed for example that they may be important in 
the dynamics of the formation of so called "ridge" in heavy ion collisions [12]. Given that the eikonal approximation 
has been extensively used in the study of hadronic scattering at high energies, it is natural to ask whether these 
longitudinal fields are accounted for in the eikonal approximation, or whether they are purely non eikonal effect. If 
the latter were the case, it would cast serious doubt on the applicability of the eikonal approximation to the nucleus- 
nucleus collision. It seems to be an implicit and unspoken assumption of the community at present, that longitudinal 
fields are indeed outside the eikonal approximation, even though this question has not been studied explicitly. We 
point out in this paper that the actual situation is in fact the opposite. By analyzing the structure of the hadronic wave 
function immediately after (eikonal) scattering we show that the longitudinal fields arise naturally and unavoidably in 
the "wake" of the projectile. The reason for this is that the Weiszacker- Williams gluons of the projectile do not gauge 
rotate "in sine" with the valence color charges in the projectile wave function. This disturbs the balance necessary 
to preserve the absence of the longitudinal fields. Thus eikonal scattering takes properly into account this important 
piece of dynamics. 

EIKONAL APPROXIMATION IN THE SYMMETRIC GAUGE. 

We are aiming to study the scattering of two hadrons, where one of them (B) is moving right and the other one 
(C) moving left with high longitudinal momentum. Most of the partons in the wave functions of those two objects 
therefore have very large positive or negative longitudinal momenta. To discuss the scattering in a symmetric manner 
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we need first off to choose a symmetric gauge. The kinematics of the process suggests that we choose 

^o=0 (1) 
In this gauge the theory is defined by the Hamiltonian density 

n=\(Tl l ? + \{F tj f (2) 

with the constraint 

D ■ n = (3) 

where II" (x) are momenta (color electric field) conjugate to the vector potential Af(x), the color magnetic field 
Ffi (x) = d t A a - djAf - gf ahc A\A c v and Df = d l 5 ah - gf ach A c t . 

We also impose a subsidiary gauge fixing condition to fix the gauge transformations that do not depend on time 

A 3 (x,i = 0) = (4) 

Imposition of the subsidiary condition in the Hamiltonian formalism is not frequently discussed. We have provided 
Appendices A and B to explain how this is done conceptually and also technically. The bottom line is that for our 
purposes, as long as we work in the classical approximation where the fields associated with the incoming hadrons are 
large, the subtleties of the procedure are not relevant, and we can treat eq.(4) as operatorial constraint which can be 
used to eliminate A 3 (x). 

Fixing A${t = 0) = 0, we can solve for II 3 from the Gauss' law: 

n 3 = -l^-n (5) 

C>3 



and the Hamiltonian becomes 

I(tf)» + i(tf«)» + i(F«) a + i 



H = o( n ') 2 + .(n 3 ) 2 + -A p ijf + o(^) 2 (6) 



where from now on we take the index i to take two values i = 1,2. Our residual gauge fixing is not complete as it 
does not fix gauge transformations which do not depend on x 3 . Thus the theory still has a constraint 



dxsDilU = (7) 

The fields and momenta satisfy the canonical commutation relations 

[^(x),nj(y)]=ia£(x-y) (8) 
and can be expressed in terms of the canonical creation and annihilation operators a* and a\ as 

/fj3 L 1 
^3 -j= K(kK k ' x + aj°(k) e - ik ' x ] (9) 



n?(x) 



; / K(k)e Jk ' x + a\ a (k)e-^] (10) 



The canonical commutation relations are as always 

[ a (k), a t ( p ) ] = (2^) 3 ^ ( 5 3 (k-p) (11) 

For simplicity of notation we will omit the color and rotational indices on the fields in the following. 
We now artificially split the vector potential A into the sum of three fields in a manner similar to [5] 

A = B + C + A (12) 
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where 

rco d 3 k 1 



B(x) = A ^^ k \- a{kykx+a]{k)e ~ lkx ^ (13) 

/—A j3l i 

C(x) = / ^ = [a(fc)e lfcx + a t (fc)e- 4fca; ] (14) 



2n y/2LUk 
~ f A d 3 k 1 

A(x) = | a ^ -j= [a(k)e- + a t (k)e --x] (15) 

where the limits in the above integrals refer to the variable k 3 . The reasoning behind this splitting is the following. 
We do not intend to discuss all states in the Hilbert space of the QCD Hamiltonian, but only those states which 
contain fast moving hadrons. The wave function of such a hadron contains overwhelmingly gluons with very large 
(positive or negative) momentum k 3 . Thus the splitting in eq.(12) corresponds to separating the degrees of freedom 
describing the right moving hadron (fields B) and the left moving hadron (fiields C). The field A contains the soft 
gluons which necessarily appear in the final state of the scattering and thus generally have to be kept. Technically 
they are also necessary to satisfy the Gauss' law eq.(7) 

The longitudinal cutoff A is arbitrary and its only role is to indicate that the longitudinal momenta in the fields B, 
C and A are vastly different. We will be thinking of it as being of the order of typical transverse momentum kj_. Since 
most of the momenta in B and C are much larger than A, we can approximate the energy u>k entering eqs.(13,14,15) 
by u)k ~ \k 3 \. We thus have 

f°° r1 3 k 1 

B(X) = A ^7m la(k)e +aHk)e 1 (16) 
c(x) = I^vm latkyk ' +a,{k)e ~' k ' ] (17) 

~ f A rfik 1 

m = | A ^- 7 ==Kk )e - +a t( k)e -- ] (18) 

Unsurprisingly, the expressions for B and C now become exactly the same as in the appropriate light cone gauges - 
the A + = for the field B and the A~ = for the field C. It is not particularly surprising, since in our gauge both 
conditions hold at time zero. Splitting the conjugate momentum in the similar way we have 



with 



n = n s + u c + u x (19) 
n B (x) = -i^°°^yB[ a(A; ) e fc_ a t (fc)e -ifcx] (20) 

IIc(x) = -if A ^-^M[a(k)e^-a^k)e-^] (21) 



00 

A 



The above equations give 



n i( x) = , I l^^f [a(k) C *--at(k) C -*-] (22) 



n = -d 3 B + d 3 c + n~ (23) 



Our next step is to express the Hamiltonian in terms of the fields B, C and A again in the approximation where the 
longitudinal momenta in the fields B and C are much greater than the transverse momenta. Also we are interested 
mainly in the case when the intensity of the fields B and C is large. More specifically we are gearing up to the 
"classical" situation where the number of gluons in each of the incoming hadrons is of order 1/g . 



f A 1 

dk 3 d 2 k ± a\k)a(k) - / dk 3 d 2 k ± a^ (k)a(k) - — (24) 
J —00 9 
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This defines our parametric counting. Assuming that the integral over longitudinal momenta is saturated on some 
large momentum scale k, we have 



a(k) ~ a\k) 



1 



Ik gk\ 



(25) 



Using this counting and keeping only the leading (order n) and the first subleading terms (order k°) in the Hamiltonian 
we find 



H^H + H 1 



with 



and 



Ho 



I 



d 6 x 



(d 3 B) 2 + (d 3 cf 



Hi = ^) 2 + -^ 3 ) 2 + - 2 {d 3 A l f + 



d. 



{DiUi + gf abc d 3 B b Bf + gf abc d 3 C b C?) 



(26) 



(27) 



(28) 



where we have dropped the tilde over A and IT; and D denotes covariant derivative in the soft field A. 
The residual Gauss' law eq.(7) becomes 



dxsDiU, + J+ + J- 







with 



J+( Xi ) = g J dx 3 f abc d 3 B b B<r; J~ = g J dx 3 f abc d 3 C^ 



(29) 



(30) 



Our goal now is to calculate the eikonal S'-matrix for the scattering process of two large hadrons (nuclei) in terms 
of the light cone currents J+ and J~. As the first step towards this goal we will rewrite the Hamiltonian in the 
interaction picture. 



THE INTERACTION PICTURE HAMILTONIAN. 



The interaction picture of quantum mechanics is defined with respect to some free Hamiltonian. In our case it is 
natural to choose as a free Hamiltonian Hq as it is the leading operator in the high energy limit. 
Any operator O in the interaction picture is defined as : 



0(t) = U*(t,t )0(to)U (t,t ) 
where Uq is the evolution operator, and H is the free Hamiltonian: 

U (t,t ) = e-^C-M Ho = [ [(dBf + (dC) 2 ] 

J X 

The free Hamiltonian Hq is just the Hamiltonian of free left and right moving energetic partons 



(31) 



(32) 



H B = I 

Jk 

H C = [ 

Jk 



Hb + He 
2?r 



2tt 



|Ae 3 | a+(As) ■ a(fe) 



\k 3 \a\k)-a{k) 



(33) 
(34) 



The interaction picture creation and annihilation operators are 

a(M) = a(fc)e~ i|fe3|t 
at(M) = a f (fc)e l|fe3|t 
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Since H does not depend on field A, the interaction picture and the Schroedinger picture A and it are identical. 
To calculate the interaction picture Hamiltonian Hi we need to look at the terms of the type 

B(-k - p,t) ■ B(p,t) (35) 

with momentum k 3 ~ k±. For the Schroedinger operators we could neglect k 3 in this expression, k 3 + p 3 w p 3 , since 
expanding in k 3 necessarily meant expanding in k 3 /p 3 . Now however there is an additional dimensional parameter in 
the game, time t. We thus have to be careful not to neglect terms of the type k 3 t which do not have to be small. 
Keeping this in mind we write 

B(-k-p,t)-B(p,t) = ^ {e{-p 3 ) a{-k-p)-a\-p) e-' l{ \ ks+ ^-\ p ^ (36) 
2|f>3| L > 

For k 3 < p 3 , 

\k 3 +Pz\ - \Pz\ - h sign(p 3 ) (37) 
Referring to eq.(28) we see that we need to calculate 

^- 9 r hc d 3 B\Bi (38) 

or in momentum space 

/1 /*oo 
dp^gf abc B\{-k 3 - P3 )Bt(p 3 ) = g-e ik ^f abc ^ dp 3 af(p)aU P ) (39) 

Analogous expression holds for the term involving the fields C with the substitution k 3 — > — k 3 . All said and done we 
obtain the "interaction" Hamiltonian in the interaction picture as 



Hl= 2 



— (A • n 4 + s(x 3 - t)j+ + 5(x 3 + t),r 



2 + n? + ^. 



{d 3 A t f 



(40) 



This expression is indeed what one naively expects. Since the most energetic "valence" partons are very fast, in the 
leading cikonal approximation their Hamiltonian is free. On the other hand the soft modes of the field A in the 
leading order are eikonally coupled to the right- and left- moving thin pancakes of color charge associated with the 
fast modes B and C. 

Our aim now is to calculate the S-matrix in the parametric regime J + <~ J~ <~ ^. In this classical regime the 
S'-matrix is given by the semiclassical expression 



5 = cxp 



(41) 



The time integral in this expression is only over the interaction time, which in the eikonal approximation is in- 
finitesimal. Thus to calculate the S'-matrix we need to find the classical solution of equations of motion only up to 
infinitesimaly short time after the interaction. The equations of motion that follow from the Hamiltonian are 



A* = 
tl a = 



1 



(42) 



gf abc Il b -2 [D-U + S(x 3 -t)J+ + 5(x 3 + t) J"] ' 

0' 



{DjFji) 



d 2 3 A« 



For t < the solution for given J+ and J is 

AL = -e(-x 3 + t)b+ a > -e(x 3 + t)b- a ' 



Ka - S{-x 3 + t)b+*-5(x 3 + t)b- a ' 



with the "classical" fields bf defined by 



dibj = J+ ; dibr = J- 
F y '(6+) =F I3 (b-) = 



(43) 
(44) 



(45) 
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Since the time integral in eq.(41) is over the infitcsimal interval, only discontinuities of the type 6{t) in A and 
singularities of the type d(t) in II can give finite contributions. It is easy to see that the only discontinuities arise 
from the extension of eq.(43) to positive times. We write the solution for arbitrary time as 



n?(x,t) = u a 0t ( x ,t) + suy(t) 

Substituting this into the equations of motion (42), we obtain 



5 A, (x,i) 
<jn°(x,t) 



51% + Df(A )^ D bc (A ) ■ <5IT - gf bcd 5A c Ii d 



i 



- g r c n D 0l ^[D ccl (A )-su a 



gf cde SA a n e 



dlSA? 



D(A )5 lj - Di(A )Dj(A ) 



5A b 



These equations have to be solved with vanishing initial conditions 

SA(t = 0) = 5U(t = 0) = 



(46) 
(47) 



(48) 
(49) 

(50) 



Examining the right hand side of eqs. (48,49) we see that they do not have singular terms. The structure of these 
equations is consistent with the solutions being regular functions. We thus conclude that SA has no 9(x 3 ± t) type 
terms and 511 has no d(x ± t) type terms. 

The classical action therefore can be calculated simply using eq.(43). The resulting S*-matrix is 



q J d 2 x b^-b i 



(51) 



This is the same expression as derived in [5] in the path integral formalism. Note that if one of the color charge 
densities is small, this expression can be simplified. Let's take for example J+ <~ g. Then to leading order in g we 
have bf = §kJ + - Substituting this into eq.(51) and integrating by parts we find 



5" = e* J d2xd2 v J+ ( x ) -ffi( x -v) J (v) 
which is the expression used in [7] for dipole-dipole scattering. 



(52) 



THE EIKONAL APPROXIMATION AND THE LONGITUDINAL FIELDS 

The point made in several papers recently[ll], is that the solutions of classical equations after collision contain 
longitudinal electric and magnetic fields. This is significantly different from the fields incoming into the collision as 
those are purely transverse. This feature is of course reproduced by our analysis of the previous section. A short time 
after the collision the solution (up to small corrections) is given be eq.(43). Calculating the longitudinal electric field 
we find 

E« = ±- [D ■ II + 5{x 3 - t)J+ + 5{x 3 + t)J~] = -gf abc b+ b (x)br c (x)9(x 3 + t)6(-x 3 + t)9(t) 
o 3 

B a 3 = \eijF* = -gr bc e^b+ b (x)bj c (x)0(x 3 + t)0(-x 3 + t)6{t) (53) 

Thus, immediately after the collision, the space between the receding fast particles is filled with the longitudinal 
electric and magnetic fields. 

Since the realization that the longitudinal fields arise after the collision, it has not been entirely clear whether 
cikonal approximation takes these effects into account. Our discussion in the previous section illustrates that this 
is not the case, as the longitudinal fields do indeed arise in the eikonal approximation. Conventionally however the 
cikonal approximation is formulated in a different frame. In this frame only one of the scattering objects (" the 
projectile" ) is fast moving, while the other one is represented by a distribution of static color fields which provide the 
target for the scattering of the partons from the projectile. It is instructive to understand how the eikonal evolution 
in this target rest frame also leads to existence of longitudinal fields shortly after collision. 

The dynamics of the collision in the target rest frame is the following. The right moving projectile starts up with 
the wave function of the form: 



\^) in =£l[a,a\j + ]\i ) ) J+ 



(54) 
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Here \ip)j+ is the wave function which depends only on the valence charge density, while the unitary operator Q 
depends on the soft gluon degrees of freedom a and and also on J + as a parameter. The soft gluon degrees of 
freedom in this approach play the role analogous to that of the soft field A of the previous section, while the charge 
density J + is due to the fast modes, the analogs of the field B. There are no dynamical modes C, as they are 
represented by the distribution of static color fields a. The operator Q is a Bogoliubov type operator as discussed 
in [10]. In order to simplify our discussion however, we will consider a classical approximation where it reduces to a 
coherent operator 



n = cxp ji J dr i d z xb a l {x)\a a l {f 1 ,x)+a { 1 a {ri,x)\ [ j (55) 

where the "classical field b as before, is related to the color charge density via 

d l b a l {x) = J+(x) (56) 

and crucially is two dimensionally a pure gauge 

dib? - djb? - gf^b^x^x) = (57) 

The operators a){rj) and 0(77) are creation and annihilation operators of soft gluons at rapidity 77. 
The operator f2 when acting on the soft gluon vacuum creates the soft fields of the form 

A1 = -b1{x)6{-x-) (58) 

The transverse electromagnetic field is: 

F+* = bl(x)S(x-) (59) 

The longitudinal electric field in the light cone quantization is not an independent degree of freedom, but is expressed 
in terms of the transverse components of the vector potential as 

E * = j+ [ D * F+l - J+S ( x ~)} (6°) 

This solves the Gauss' law in the light cone gauge. With the transverse field of eq.(59), the eq.(56) ensures that the 
longitudinal electric field before the collision vanishes. The longitudinal magnetic field obviously vanishes as well since 
bi is two dimensionally a pure gauge. Thus the incoming projectile contains only transverse electromagnetic fields 
which are concentrated within a "shock wave" - the transverse plain moving along with the valence charge density 
created by the fast partons. Those are of course the well known Weizsacker - Williams fields. 

The target in this approach is represented by an ensemble of color fields a a . As discussed at length in [4], [9], [3] the 
projectile emerges from the interaction region with the wave function 

\^) out = n[Sa,SalSJ+}\^) S j + (61) 
where the single gluon 5-matrix S is expressed in terms of the target field 

S*(x)=cxp{i J dx-T a a a (x,x- )} (62) 

For the operator O of eq.(55) this gives 

rt[Sa,Sa\SJ + ] = expji J drjd 2 xb? (x) a? (17, x) + a\ a (r], x) J (63) 

with 

b? = S^[SJ + ] (64) 
The transformed operator f2 now creates the transverse field 



F+* =bfa)6(x-) 



(65) 
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This at first sight looks very similar to the fields before the collision. There is however a crucial difference. These 
transverse fields accompany the rotated color charge density SJ + . Obviously the color charge density is not rotated 
entirely "in sine" with the field in the sense that 

dh ^ SJ + (66) 

Therefore Gauss' law requires nonvanishing of the longitudinal electric field. The longitudinal electric field after the 
scattering is indeed given by 

E 3 = ^ T [D l F+ l -SJ+5(x-)] = 9(t-x 3 )d l ^b l [SJ + }-b l [SJ + ]) (67) 

and is clearly nonzero. Likewise, since 6, is not a two dimensional pure gauge, the longitudinal magnetic field does 
not vanish cither: 



d t b a 3 - djb? - gf aM b1{x)b c 3 {x) 9{t ~ x 3 ) (68) 



One note of caution is that eqs. (67,68) are valid only close to the longitudinal position of the projectile - that is only 
between the projectile and the target. When extended beyond the longitudinal coordinate of the target, the integral 
1/9+ will pick up additional contributions due to the target fields, so that presumably the longitudinal fields vanish 
for negative x 3 . Otherwise eqs. (67,68) closely parallel the structure of longitudinal fields in the symmetric frame. 

APPENDIX A. GAUGE FIXING IN THE HAMILTONIAN FORMALISM. 

In this Appendix we discuss how to consistently impose a gauge fixing for the residual gauge symmetry in the 
Hamiltonian formalism. Consider a system with a Hamiltonian H which is invariant under the action of some 
symmetry generators C a : 

[H,C a ]=0 (69) 
Also the space of physical states is restricted to satisfy 

C a \*/J)=0 (70) 
All observables in such a theory are necessarily also gauge invariant 

[C a ,O}=0 (71) 
Thus an expectation value of any observable in any physical state is given by the formal expression 

(O) = J dxidy a H>*(x,y)0*(x,y) (72) 

Here we have intentionally separated the coordinates into two sets {xi} and {y a }, since we intend to impose gauge 
fixing conditions and solve them for y a . We want to impose the gauge fixing 

G a (x,y)=Q (73) 

Let us a la Fadeev-Popov multiply the matrix element in eq.(72) by unity 

d\ a S[G a (x x ,y x )]\detj^\ (74) 

where x x and y x as usual are x and y transformed by the gauge transformation with gauge parameters A. Just like in 
the usual Fadeev-Popov approach we can now change variables in the integral x, y — ► x x , y x . Due to gauge invariance 
of both the wave function 'J and the operator O, the dependence on A drops out of the integrand, and the A integral 
gives an irrelevant constant. We thus have 

/SC 
dx l dy a ^*(x,y)0^(x,y)S[G a (x,y)}dctj^\ 

= J dx^*{x,y = Y{x)) [OV(x,y)) y=Y{x) \det~ lS -^det 6 -^\ y=Y(x) (75) 
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Here y a = Y a (x) are solution of the gauge fixing conditions eq.(73). This equation has a simple interpretation. Any 
matrix element of the type of eq.(72) can be calculated as a matrix element on the smaller Hilbert space, the space 
spanned by coordinates X{ only 



On this space the norm and the scalar product are defined with a nontrivial measure 



fj,(x) 



Idct 



-l6G 



SG n 



5yp 5 A/3 



det 



\V=Y(x) 



(76) 



(77) 



The wave functions and the operators are related to those of the original gauge invariant theory. For the wave 
functions this correspondence is simple to define 



\f/(x) = if>(x,y = Y(x)) 



(78) 



For the operators the correspondence is more subtle. It is straightforward if the operator O depends only on coordinates 
and not on the momenta. In this case, just like for the wave function 



0(x)=O(x,y = Y(x)) 



In general the correspondence is given by 



0(x,pWx) = [0*(x,y)] y=Y(x) 



(79) 



(80) 



The nontrivial statement is that such a correspondence exists for any operator O so that the operator O is the same 
for all gauge invariant wave functions , 3>. 

In general finidng this correspondence requires a nontrivial amount of work. We will not delve into the discussion 
of the most general case, but instead will only consider a simple situation when the constraints C a as well as gauge 
conditions G a are linear functionals of the coordinates and momenta. In this case the constraints C a can be solved 
simply for the momenta conjugate to y a 



Pa 



9a{Xi,yf3,Pi 



(81) 



We first show how Pi - momenta conjugate to x i: are represented when they act on <3>. To do this we have to remember 
that iff is only a function of combinations of coordinates Ri which commute with C a . We can then write 



dV{x,y) _ d$> dRj 



dxi dRj dxi ' dxi dRj 

On the other hand the gauge invariance of R means 



\y=Y(x) 



dRj 

dxi 



9* 

or~ 



y=Y(x) 



dRj dRj dY a 
dxi dy a dxi 



[Ra,P/3 ~ gp(xi,yp,Pi)} = i 
Combining these equations we find 



d^(x,y) , 
dx t 



\y=Y(x) 



( dR a dR a dgp 
\dy f3 dxj pj 

Mi .dm 
13 dxj 



= 



ith 



dY a dg a 
dxi dpj 



(82) 



(83) 



(84) 



(85) 



By the generalization of this argument we can show that the same is true for the second derivative of iff provided the 
matrix M is a c- number matrix (which is where the linearity of C and G is important). 



d 2 if{x,y) i 
dxidxj 



\y=Y(x) 



M lk M,i 



d 2 ^{x) 
dxkdxi 



(86) 
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Thus we see that on the reduced Hilbert space the commutators of the coordinates Xi and momenta pi are modified 

[xi,pj] = -iMij (87) 

The matrix M is in fact the Dirac bracket as shown in Appendix B. 

The momentum p a can be eliminated in favor of pi using the constraint eq.(81). One has to be careful though, as 
this constrained can only be used when the operator p a acts on a gauge invariant state. Thus 



Pa\^) = 9a(x,y,Pi)\ y =Y(x)\^)- 

1 

2 



(88) 



PaPpl^) = Pa9/3{x,y,Pi)\^) 

P a Pi\^) = Ptg a (x,y,Pt)\y) 



{gp(x,y,Pi),g a (x,y,Pi)} - i (j^p + 

y=Y {x) \y) 



\y=Y(x) 



|M>) 



g a {x,y,Pi)Pi - i-^r— 

OXi 



Thus we get the following algorithm to represent an arbitrary operator O on a gauge fixed Hilbert space.. 

1. Order factors of momenta and coordinates such that all the momenta are to the right, and all the coordinates 
to the left. 

2. Carefully realize the Gauss law constraints using eq.(88) and its generalizations if necessary. 

3. Represent all momenta pi by derivatives of In the classical approximation this amounts to replacing the 
Poisson brackets of pi by Dirac brackets. In the quantum theory this requires more work, in particular new quantum 
terms arise which involve derivatives of the Dirac brackets. When constraints are linear, the Dirac bracket procedure 
is exact also in the quantum case. 

4. Finally express all y a in terms of X{. 

This concludes our discussion of gauge fixing in the Hamiltonian formalism. Finally we note that although the 
Gauss' law in QCD is nonlinear, for the purpose of the present paper the discussion in this appendix suffices. The 
reason is that our aim was to derive the expression for the ^-matrix for large fields bf . In this situation it is sufficient 
to expand the Gauss' law to first order in fluctuations around the classical fields. In leading order the constraint 
becomes linear and the procedure described in this appendix applies. Moreover, imposing Dirac brackets is equivalent 
to operatorially solving the gauge fixing condition, which justifies the quantization procedure used in this paper. 
Beyond the (semi)classical approximation the qunatization procedure is well defined, but is more subtle as explained 
above. 



APPENDIX B. EQUIVALENCE WITH DIRAC BRACKETS. 

Here we show that the commutation relations discussed in Appendix A are equivalent to Dirac bracket. Consider 
the Dirac bracket quantization of the system with constraints 

C a = p a - g a (xi,yi,pi) = 0; G a = y a - Y a (x) = (89) 

The Dirac brackets are defined as 

{A, B} D = {A, B} P - {A, c a } P C-^{c b , B} P (90) 
where c a is the complete set of second class constraints (in our case C a and G a ) and the matrix C a b is defined as 

Cab = {Cc,C h } P (91) 

The matrix C is antisymmetric. In our case it has vanishing diagonal blocks, while the off diagonal blocks are 

dg a dYp 



L a ,p = {c a , g (3 }p = (s afi + I^H) 



dpk dx k 

Now consider the Dirac brackets between the coordinates and momenta 

dg a dYp dg a dY 1 dg 1 dYp 



f , x dY a !%j , . dY c 



dx. 



dp k dx k + dp k dx k dpi dx t 



dgp 
d Pj 



(92) 



(93) 
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This has to be compared with the result obtained in the text 



{Pi > j 



-iMij 



i\S - - d9a ] 
dxi dpj 



g _ dYadga 9Y a dg a dY dgp 
13 dxi dpj dxi dpk dxk dpj 



(94) 



The two expressions obviously coincide. We thus see that our analysis in Appendix A reproduces the Dirac bracket 
quantization procedure. We note however, that the equivalence with the Dirac bracket quantization only holds on the 
classical level, or in a fully quantum theory when all the constraints are linear. If the constraints are nonlinear, then 
the quantum theory involves the Faddeev-Popov determinant and the additional contribution to the measure, which 
can not be obtained classically. Also if the matrix M depends on fields (which is the case for nonlinear constraints) 
one has to be careful with the operator ordering and order operators consistently following the technique discussed 
in Appendix A. 
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